Dynamics of fullerene coalescence 
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Fullerene coalescence experimentally found in fullerene- embedded single-wall nanotubes under 
electron-beam irradiation or heat treatment is simulated by minimizing the classical action for 
many atom systems. The dynamical trajectory for forming a (5,5) C120 nanocapsule from two C60 
fullerene molecules consists of thermal motions around potential basins and ten successive Stone- 
Wales-type bond rotations after the initial cage-opening process for which energy cost is about 8 
eV. Dynamical paths for forming large-diameter nanocapsules with (10,0), (6,6), and (12,0) chiral 
indexes have more bond rotations than 25 with the transition barriers in a range of 10-12 eV. 

PACS numbers: 61.46.+w,34.10.+x,81.07.-b 



Carbon nanostructures spanningfrom fullerenes 0] to 
single and multi-wall nanotubes [3, nano-peapods Q, 
and "Y" and "T" junctions Q have attracted much at- 
tention from scientists for last two decades due to their 
unique structural, mechanical, and electrical properties 
0. In particular, carbon nanotubes and related junc- 
tions have great potential for applications on nanometric 
electronic devices such as quantum functional transistors 
and rectifiers since nanotubes can be semiconductors or 
metals depending upon their chiralities and diameters . 
At the moment, however, development of the nanotube- 
based technology is limited because delicate control of 
their structures is not yet feasible. 

Recently, merging processes between carbon nanos- 
tructures have been proposed as a new synthesis tech- 
nique to construct more complicate functional structures: 

(1) Two nanotubes coalesce into diameter-doubled one 
in a high-energy (1.25 MeV) electron beam condition 0. 

(2) In the same condition, initially crossed nanotubes are 
found to form some functional junctions shaped like "Y" 
and "T" characters Q. Since each carbon atom is dis- 
placed approximately every 100 sec by knock-on colli- 
sions with the high-energetic electrons, the irradiation- 
induced vacancies may play a crucial role in these merg- 
ing processes. Terrones and his coworkers 0, studied 
the vacancy-related merging mechanism utilizing the or- 
dinary molecular dynamics technique. Though introduc- 
ing various functional structures, this high-energy pro- 
cess may be not appropriate for nanometric device ap- 
plication because of the uncontrollable defect formation. 

(3) A low energy growth technique for nanotube bun- 
dles, where the chiral indexes of all individual nanotubes 
are almost perfectly mono-dispersed, was proposed in the 
mixture of Cqo molecules and catalyst Ni pillars with 
a help of a ~1.5 T magnetic field |8j. The annealing 
temperature is just 950 °C. (4) In another low energy 
situation, a single-wall nanotube can be synthesized by 
merging fullerene molecules inside a nanotube vessel in 
nano-peapod systems 0, 0] . Of special importance is 



the last technique of fullerene coalescence because it is 
a catalyst-free and bottom-up process, having a poten- 
tiality for diameter and chirality control of single-wall 
nanotubes. 

Nano-peapods Q consist of array of fullerene molecules 
(inner peas) and a single- wall carbon nanotube (an outer 
pod), and all components are separated from the others 
by the van der Waals distance. The fullerene molecules 
as an one-dimensional crystal are spaced by about 0.97 
nm or less, and the filling fraction is over 90 %. The high 
fraction is associated with the exothermic encapsulating 
process 0,^1- Depending on preparation methods, the 
various fullerene species can be Cgo, C70, Cs2, and met- 
allofullerenes 0, U E3 

Electron-beam irradiation Q or heat treatment 0] 
on these carbon composites induces coalescence of core 
fullerenes into an inside single- wall nanotube, of which di- 
ameter is 0.71 nm smaller than that of the outer 'vessel' 
tube. Since the electron beam in this case is irradiated for 
hundreds of seconds and its energy is about 0.1 MeV, we 
can not expect the irradiation-induced vacancies as like 
Refs.0and[3 Actually, the outer wall was little damaged 
during the coalescence process by the electron beams |9j . 
Iijima and his coworkers baked nano-peapods up to 1200 
°C for 14 hours. They found that the inside fullerenes 
start to coalesce at 800 °C and complete transforming to 
a single-wall nanotube at 1200 °C. They also proposed 
that diameters of the synthesized nanotubes are varying 
from 0.6 to 0.9 nm depending on the size of the outer 
'vessel' tubes in a range of 1.3-1.6 nm 0. This low en- 
ergy merging process of fullerenes have great implication 
on diameter and chirality control in nanotube growths 
as we mentioned above, but little was known about its 
dynamical process. Furthermore, the ordinary molecular 
dynamics study for the process is not applicable because 
it never happens within the simulation time in use of the 
modern supercomputers. 

Instead, topological consideration was made by Ueno 
et al. to the case of a C120 capsule with (5,5) chi- 
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ral index from two dimerized Cgo fullerenes |14j . They 
found 22 successive Stone- Wales-type bond rotations as 
a topologically acceptable pathway for the transforma- 
tion. Recently, hypothetic transition states between two 
consecutive local minima of the pathway was examined 
(lif. Similar topological analysis for head-to-head coa- 
lescence of capped nanotubes and fullerenes was made 
[13: (5,5)+(5,5), (10,10)+(10,10), and C 60 +(10,10). 
However, these do not guarantee real dynamics of the 
coalescing process beyond topology. Fullerene colli- 
sion 0, and nanotube collision 0] were dealt 
with molecular dynamics techniques with the assump- 
tion of one constituent's high energetic incidence on the 
other. Fullerenes packed well inside single- wall nan- 
otubes (nano-peapods) are not relevant with such a high 
energetic incidence. 

In this Letter, to simulate the low-energy merging pro- 
cess of fullerene molecules with appropriate molecular dy- 
namics information, we use the action-derived molecular 
dynamics method, recently proposed by Passerone and 
Parrinello l2(l| with some modifications for many atom 
systems [21|. 

The discretized version [2^ of the classical action (S) 
for numerical simulations is given by 

p-i 

S=J2&(Kj-Vj), (1) 

3=0 

where P is the total number of steps for discretization 
of simulation time r, A = t/P is the constant time in- 
terval, and Kj (Vj) are the kinetic (potential) energy of 
the configuration at the j-th time step. The instanta- 
neous velocity is approximated as Vj = (q^ — q J+ i)/A, 
and then Kj = mv| /2 where m is the mass of each par- 
ticle. Dynamical paths can be obtained by finding the 
stationary point of the action, i.e., SS=0, but in general 
it is a very difficult task to find the stationary point be- 
cause it involves a root-finding problem which requires 
the second derivative of the potential. Furthermore, the 
action is not bounded and the stationary point can be 
either a local minimum/maximum, or a saddle. 

To obtain effectively a dynamical path that is globally 
close to the Verlet trajectory, i.e., an optimal solution for 
the discretized action, we employ a rather simple min- 
imization technique, where the action is supplemented 
with some penalty functions ; 

p-l 

Q(qj,E; K) = S + fi E J2 ( E i ~ E ? 

3=0 

P-l P-l 

+ MP J2 p/ + ml E V + »k «*i> - K f • ( 2 ) 

3=0 3=0 

Here, P., (ly) is the instantaneous total linear (angular) 
momentum, and Ej = Kj + Vj is the instantaneous to- 
tal energy of N particles at step index j. E and K is 




FIG. 1: (color online) Initial and final states for forming var- 
ious C120 nanocapsules as precursors for fullerene-based nan- 
otube growth. 

the target total and kinetic energy, respectively, and all 
the /x's are numerical parameters that enforce the corre- 
sponding constraints. The (• • •) means the time average. 
With he, A*p, and /iL, we enforce the conservation laws 
of total energy, and total linear/ angular momenta on the 
optimized trajectory. The last term in Eq. (2) is asso- 
ciated with the virial theorem that classical trajectories 
should satisfy, causing thermal motions around potential 
basins; in thermodynamic language, the virial theorem 
corresponds to the equipartition theorem. 

In our action-derived molecular dynamics simulations, 
we minimize the generalized action O with respect to 
various atomic paths using a conjugate gradient tech- 
nique for the given initial and final iV-atom configura- 
tions, as shown in Fig. 1. The potential energy and 
Hellmann-Feynman forces are calculated using a tight- 
binding model for carbon systems Details of formu- 
lation and numerical procedure will be presented else- 
where mi. 

For the fullerene coalescing process, we can reasonably 
assume that its initial state is two Ceo molecules sepa- 
rated by ~1 nm However, for the final configura- 
tion, we can expect various C120 nanocapsules with differ- 
ent diameters and chiralities because experimentally the 
diameters of the synthesized inner tubes ranges in 0.6-0.9 
nm The simplest case is the C120 nanocapsule with 
diameter of 0.7 nm and chiral index of (5,5) [Hereafter 
we call it the (5,5) nanocapsule], as other workers con- 
sidered 0, flol flq | . We also present possible diameter 
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FIG. 2: Potential energy variation along the obtained dynam- 
ical trajectories of fullerene coalescence for forming various 
C120 nanocapsules. In (a) and (b), the solid (dotted) lines 
indicate the potential (total) energy of the (5,5) and (10,0) 
capsules, respectively, along the time step. We divide the tra- 
jectories into three parts as reactant, activation, and product 
regions in series. In (a), ten small arrows correspond to ten 
Stone- Wales- type bond rotations. In (b), the single arrow in- 
dicates when the intermediate (5,5) capsule arises. In (c), we 
depict the potential energy variation in case of the (6,6) and 
(12,0) C120 capsules. Also, the small arrow indicates forming 
the intermediate (6,6) capsule in the trajectory of the (12,0) 
capsule. 



thickening process by forming (10,0), (6,6), and (12,0) 
nanocapsules with diameters of 0.8-9.8 nm, as shown in 
Fig. 1. Assuming that all the atoms move least dis- 
tance from initial to final states in order to accelerate 
the numerical convergence, we choose the initial orienta- 
tion of two fullerene molecules differently depending on 
the cap structure of the final nanocapsules. With this we 
can remove free rotational motion of a single molecule 
in the full molecular dynamics simulation, with saving 
much computation time. 

We perform action-derived molecular dynamics simula- 
tions using the parameters of /Lt£=10 8 , ^,- p =H'l=^k =10 , 
P=199, and t=P A— 1.2 ps, and display the results in 
Fig. 2. With the numerical parameters, we find that 
the total energy is well conserved over all time steps for 
all the nanocapsules, whereas the linear and angular mo- 
menta is well converged for the (5,5) and (10,0) capsules, 



but less strictly for (6,6) and (12,0) capsules. Since E 
and K cannot be independent variable in a microcanon- 
ical problem where total energy of systems is conserved, 
we relax K within the minimization technique, regarding 
it as an extra variable after achieving some convergence. 

For the coalescing process of the two Ceo molecules into 
the (5,5) C120 nanocapsule, we can divide the obtained 
dynamical trajectory [See the supplementary movie] into 
three parts as reactant, activation, and product regions 
in series, as in Fig. 2(a). In the reactant region around 
the initial configuration, the 120 atoms of the two Ceo 
molecules are randomly vibrated around their energy 
minima with preserving fullerene cage structures. This 
motion can be regarded as thermal motion, and be con- 
sistent with random trajectories around potential min- 
ima of the Mueller potential in Ref. H3- The relation 
(Kj) w (AVj), the virial theorem for the harmonic poten- 
tial, is satisfied. In the product region near the final state, 
we similarly find a thermal vibrating motion of the (5,5) 
C120 nanocapsule. The coalescence occurs in-between, 
i.e., the activation region of about 0.3 ps. At the tran- 
sition state, two fullerenes start to merge with breaking 
some C-C bonds to open cage, and then a peanut-shaped 
intermediate structure appears at step index 87. Then, 
we can see 10 successive Stone- Wales- type bond rotations 
forming small potential peaks due to the occurrence of 
undercoordinated carbon atoms. The ten rotations are 
attributed to the 5-fold symmetry of the initial and final 
states, as discussed in Ref. El The peak height is about 
2-5 eV. For the whole coalescence process, the energy 
barrier is estimated to 8.3 cV. This is quite comparable 
to the energy barrier for single Stone- Wales transforma- 
tion in a Ceo molecule |24|. 

For the (10,0) nanocapsule, we also find the reac- 
tant and product regions of the dynamical trajectory 
with randomly vibrating atomic motion around both 
the potential basins. We count about 25 Stone- Wales- 
type bond rotations in-between activation region, lasting 
about 0.6 ps to complete the transformation. An inter- 
mediate (5,5) C120 nanocapsule is found at step index 
90 on the trajectory for the (10,0) capsule, indicated by 
the arrow in Fig. 2(b). Thus, we find a transformation 
route: two Ceo fullerenes — * (5,5) C120 nanocapsule — ► 
(10,0) C120 nanocapsule. In the former part, a transi- 
tion state occurs with an energy barrier of 9.9 eV, taking 
about 10 bond rotations, and in the latter part the (5,5) 
nanocapsule transforms into (10,0) nanocapsule through 
about 15 bond rotations, with an energy barrier of less 
than 8 eV. Even though the former process is concep- 
tually same to the coalescence process for forming (5,5) 
nanocapsule of Fig. 2(a), the energy barrier is higher 
by 1.6 eV because of difference in the initial molecular 
orientation as shown in Figs. 1(a) and (b). 

In the cases of (6,6) and (12,0) nanocapsules, the whole 
processes require more bond rotations than 30. With the 
simulation time of 1.2 ps, the obtained dynamical trajec- 
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TABLE I: Fullerene coalescence into various C120 nanocap- 
sules. Diameter of capsules (D), formation energy (Hf), over- 
all potential energy barriers (Vb), average kinetic energy in 
Kelvin ({Kj)), and O values are compared for the used target 
energy (E). 



D (A) H f (eV) V b (eV) (Kj)(K) Q(a.u.) E (a.u.) 



(5,5) 


6.85 


-14.53 


8.3 


556 


599.5 


-29.85 


(10,0) 


8.02 


-15.81 


9.9 


360 


368.9 


-29.85 


(6,6) 


8.32 


-16.25 


11.4 


239 


227.7 


-29.80 


(12,0) 


9.78 


-16.30 


12.4 


353 


300.7 


-29.78 



tories shown in Fig. 2(c) entirely belong to the activa- 
tion process, different from those for the (5,5) and (10,0) 
nanocapsules in Fig. 2(a) and (b). If we use sufficiently 
large r to cover 30-60 bond rotations, we may get such 
thermal vibrating modes near the initial and final states, 
but this demands much computational cost. With the 
limited simulation time, even though we can not guar- 
anteed the real dynamics [In fact, total linear/angular 
momenta are poorly conserved], we can estimate sat- 
isfactorily activation energy barriers within the action- 
derived molecular dynamics formalism. The calculated 
values are 11.4 (step index 97) and 12.4 eV (step index 
50) for the (6,6) and (12,0) nanocapsules, respectively, 
rather high than those for the (5,5) and (10,0) nanocap- 
sules. The limited simulation time causes the simultane- 
ous Stone- Wales-type bond rotations, resulting the high 
energy barriers. Another reason is that, since we pre- 
serve the 6-fold symmetry in the initial and final states 
of these systems, double bonds of the hexagons should be 
broken to coalesce initially. This initial opening process 
needs an energy cost of about 10 eV, which corresponds 
to step index 0-30. We also find an intermediate (6,6) 
nanocapsule at step index 130 in the path of forming the 
(12,0) nanocapsule, indicated by the arrow in Fig. 2(c). 

We summarize all results and the target energy (E) 
used for these four nanocapsules in Table I. As the di- 
ameters increase, the capsules are more stable with re- 
laxing strain energy, but estimated energy barriers are 
higher. The large energy barriers over 10 eV in the cases 
of the (6,6) and (12,0) nanocapsules are attributed to 
the limited simulation time and thus simultaneous bond 
rotations. Our estimated values for transition barriers 
are quite reasonable for low-energy merging process of 
fullerenes inside single-wall nanotubes. The value of 10 
eV corresponds to system temperature of 1450 K for two 
Ceo molecules, and the experimental temperature for the 
merging process is 1500 K We also display the aver- 
age kinetic energy (Kj) used for ensuring the virial theo- 
rem, and the discretized Onsager-Machlup actions O [2^| . 
The obtained O values are rather large, especially for the 
cases of the (5,5) nanocapsule due to the thermal motions 
in dynamical trajectory. We note that the action-derived 
molecular dynamics cannot offer a unique solution for the 



dynamical paths of a specific rare event. We may obtain 
different dynamical paths with various O and (Kj), de- 
pending on the target energy, initial seed, and minimiza- 
tion schedule for the functional Eq. (2). Nevertheless, 
we believe that the 8 trajectory generally reflects the 
real dynamics of the system under consideration via the 
constraints for the classical laws. By imposing the appro- 
priate dynamical laws and virial theorem, we can obtain 
an estimate, not exact value, for transition energy barrier 
and transformation process for specific rare events based 
on the real dynamics. 

In conclusion, we investigate the molecular transfor- 
mation of two Cgo fullerenes into various C120 nanocap- 
sules, and find dynamical paths for the fullerene coales- 
cence using the action-derived molecular dynamics tech- 
nique. In the dynamical process of forming the various 
C120 nanocapsule with the chiral indexes of (5,5), (10,0), 
(6,6), and (12,0), we find that the process consists of suc- 
cessive Stone- Wales-type bond rotations with the transi- 
tion energy barrier of about 8-12 eV. Further growth 
via coalescence between nanocapsules and Ceo fullerenes 
may be natural extension of the above coalescence pro- 
cess with great similarity, requiring future study. 
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